Noise spectra of an interacting quantum dot 
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I. INTRODUCTION 

The fluctuations of the electric current flowing through 
a mesoscopic system, connected to several electronic 
reservoirs (e.g., a constriction between two reservoirs) are 
described by the noise spectrum. This spectrum contains 
valuable information on the characteristics of the sys- 
tem, which cannot be detected by dc conductance mea- 
surements. Examples include information on Coulomb 
interactions, on the quantum statistics of charge carri- 
ers, on the number of available transport channels and 
on electrons entanglement Another example comes 
from heat transport: recently it was shown that at fi- 
nite frequencies, the noise associated with correlations of 
energy currents in a mesoscopic constriction is not re- 
lated to the thermal conductance of the constriction by 
the fluctuation-dissipation theorem, but contains addi- 
tional information on the electric conductance, especially 
at very low temperatures where the heat conductance is 
vanishingly small. 5 

The non- symmetrized noise spectrum is given by the 
Fourier transform of the current-current correlation 2 - 

/oo 
dte iut (6i a (t)6I a ,(0)) , (1) 
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where a and a' mark the leads carrying the current be- 
tween the electronic reservoirs and the system at hand. 
In Eq. (Q]), SI a = I a — (I a ), where I a is the operator 
of the current emerging from the ath reservoir, and the 
average (denoted by (•••)) is taken over the reservoirs' 
states. Mesoscopic devices typically operate at low tem- 
peratures T, allowing measurements at relatively high 
frequencies, which obey fku > k B T. In this range one 
may observe quantum deviations from the more familiar 
thermal fluctuations £ The non-symmetrized noise spec- 
trum can be interpreted as due to absorption of energy 
(positive frequencies) or emission of energy (negative fre- 
quencies) by the noise source^— This interpretation does 



not hold for the symmetrized version of the noise spec- 
trum, given by C£?M = [C aa ,{u) + C aa ,(-u>)}/2. 

Measuring the non-symmetrized noise spectrum is a 
challenging task. It involves a "quantum detector", 
where the measured frequency matches an energy gap 
in the detector. Deblock et al£ used a superconductor- 
insulator-superconductor (SIS) junction as a noise detec- 
tor, and the measured frequency corresponded to the SIS 
energy gap. Later on, Gustavsson et al^ used a two- 
level system (TLS) as a detector, where the measured 
frequency corresponded to the energy difference between 
the two levels. A similar TLS system was used in Ref. [ll| 
where the equilibrium non-symmetrized noise spectrum 
was detected. More information about quantum noise 
measurements can be found in a recent review by Clerk 
et al& 

In previous work we have used the scattering formal- 
ism to show that the noise spectra of a non-interacting 
quantum dot exhibit features such as steps and dipS) 12 i 13 
whose relative locations reflect the resonance levels of the 
quantum dot and the energy spacings in-between them. 
A similar dip structure was found by Dong et alJ^ for 
two non-interacting or highly interacting dots placed on 
the branches of an Aharonov-Bohm interferometer, using 
a quantum master equation. The dip was supplanted by 
a peak once the repulsive interaction was strong enough. 
These authors also found that the depth of this dip os- 
cillates with the Aharonov-Bohm fluxj 14 confirming the 
results of Ref. [TH, which pointed out that the dip is due 
to interference between resonance levels. Later on, Wu 
and Timmi£ used the Floquet master equation to study 
the noise of a single-level quantum dot driven by either 
an ac gate voltage or by a rotating magnetic field. In 
both cases they found a peak or a dip, depending on the 
nature of the ac driving field, at a frequency matching 
the energy difference between pairs of Floquet channels 
(associated with transport channels). Recently, Marcos 
et a/*i£ used a quantum master equation to study a sys- 
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tern similar to the one analyzed in Refs. [H and[l3|, and 
found an interference-induced dip in the noise spectrum 
of a quantum dot coupled symmetrically to two leads. 
This result contradicts that of Rcf. which found that 
such a dip appears only when the leads are coupled asym- 
metrically to the dot. We explore this contradiction in 
Sec. IIV1 where we confirm that the dip disappears for 
a symmetric coupling of the dot to the leads. Moreover, 
Marcos et alM argue that the dip is immune to the appli- 
cation of a finite bias voltage and to finite temperatures. 
In Sec. IIVI we find, in agreement with the results of Ref. 
Htl that the dip location is unchanged by the application 
of a finite bias and\or temperature; its depth though is 
affected, and it may even turn into a peak. Very recently, 
Hammer and Belzigi^ extended the model of Ref. [l3| by 
considering a periodic ac bias voltage and found a similar 
step- like structure for the noise as the ones found in Refs. 
[I^ll3l and to the one found here. 

The noise spectrum of a quantum dot, as a function of 
the frequency, is indeed rather sensitive to the coupling 
of the dot to the leads J^r— This sensitivity has been 
analyzed by Thielmann et alr^ and measured by Gus- 
tavsson et al— who have monitored the shot noise as 
a function of the asymmetry of the dot-leads couplings. 
This asymmetry turned out to be in particular impor- 
tant in the experiment of Konemann et a/, 20 where the 
current through a highly asymmetrical quantum dot was 
measured as a function of the applied bias. These authors 
attributed the observed dependence of the broadening of 
the slope of the steps in the I — V characteristics on 
the polarity of the bias voltage to electronic interactions 
taking place on the dot: They argue that the bias volt- 
age determines whether the energy level, located between 
the two chemical potential of the leads, is predominantly 
empty or occupied.— We expound on this point in the 
Appendix. 

In the present paper we study the noise spectrum of the 
current passing through an interacting spinless quantum 
dot, coupled to two electronic reservoirs via two leads, 
and thus extend our previous calculations to include the 
Coulomb interaction. That interaction is treated within 
the Hartree-Fock (HF) self-consistent approximation; by 
this, the Hamiltonian is turned effectively into a single- 
particle one, enabling us to employ the scattering for- 
malism in our calculations. Those are augmented by 
a self-consistent determination of the Hamiltonian ef- 
fective energies and hopping terms among the levels. 
Having mapped the Hamiltonian to an effective single- 
particle Hamiltonian we have constricted the effect of the 
Coulomb interaction. For example, the Kondo effect, a 
many-body phenomenon associated with the screening of 
the spin of the quantum dot (for a dot with an odd num- 
ber of electrons) by the conduction electrons of the leads 
at low temperatures, can not be observed in our formal- 
ism. Zarchin et ali^ have measured the low-frequency 
noise spectrum to extract the backscattered charge from 
a two terminal system with a quantum dot in the Kondo 
regime and found it to be highly non trivial. Recently, 



using a real-time functional renormalization group ap- 
proach, Moca et al. 22 found two dips in the noise spec- 
trum that originate from the Kondo effect. 

Self-consistent Hartree and Hartree-Fock approxima- 
tions have been widely used to study the spectra of quan- 
tum dots.— One expects the HF approximation to be bet- 
ter suited for systems having large numbers of electrons 
and states^ Several authors used the HF approximation 
to calculate the occupation of levels on dots coupled to 
reservoirs. Sindel et at— considered a two-level spinless 
quantum dot, and found agreement between the results 
of the Hartree approximation and those of the numeri- 
cal renormalization-group for levels whose width is small 
compared to their separation. When that width is com- 
parable to (or larger than) the level spacing the charging 
of a level has been found to oscillate with the gate voltage. 
Goldstein and Berkovits^ used the HF approximation to 
study a multi-level quantum dot. They also obtained a 
non-monotonic behavior of the level charging when one 
level was much broader than the others; but when the 
width of the levels was small compared to their spacing 
the level charging was found to be monotonic as a func- 
tion of an applied gate voltage. Below we consider many 
resonant levels, whose width is small compared to their 
spacing. We expect the HF approximation to be valid 
in that regime. Indeed, our results for the level occupa- 
tions are in agreement with the results of Refs. and 
[26j . In the opposite regime, where the width of the levels 
is much larger than the level spacing, the so-called open 
dot regime, solving a kinetic equation for the occupancy 
on the dot, Catelani and Vavilov^ have found correc- 
tions to the shot noise due to electronic interaction, which 
do not manifest themselves in the conductance, that de- 
pend non-trivially on bias, temperature, and interaction 
strength. 

The linear-response ac conductance of a multi-level 
quantum dot coupled to a single lead was measured 
by Gabellli et al.j 2 ^ and found to exhibit a quantized 
behavior which suggested the use of this device as an 
on-demand coherent single-electron source. Expanding 
the ac conductance at low frequencies, and using the 
Hartree 2 ^ and Hartree-Fock2S approximations, the oscil- 
lations in the conductance were calculated and found to 
be in agreement with the experimental data. Here we 
extend these calculations to the case in which the multi- 
level dot is coupled to two leads. We find that the dc con- 
ductance peaks broaden as the Fermi energy is elevated 
and that the energy difference between two consecutive 
peaks grows linearly as a function of the Coulomb inter- 
action strength. 

We begin in Sec. |TT] by defining the different possi- 
ble noise spectra, and continue with a brief description 
of the scattering formalism. We then explain in detail 
the model Hamiltonian and its mapping onto an effec- 
tive Hamiltonian using a mean-field approximation. The 
results of the self-consistent approximation and the noise 
spectra are presented in Sec. IIII1 and summarized and 
discussed in Sec. IIVI In view of the controversy men- 
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tioned above, we devote a great part of Sec. [TV] to ana- 
lyze the dip appearing in the noise spectrum of a (non- 
interacting) two-level dot, and study in particular its 
modifications under the application of a bias voltage and 
finite (although low) temperatures. In the Appendix we 
explore the possible experimental detection of the Fock 
terms. 



where 

A/=(4+4)/2. (9) 

As opposed to the net current operator Eq. ([3]) , for which 
(I) in the presence of a finite bias voltage, (A/)=0. 31 
It was shown in Ref. [13| that at zero frequency the noise 
associated with AI is always zero, i.e., C'( + \Q) = 0. 



II. DETAILS OF THE CALCULATION 

A. The noise spectra of a two-terminal system 

Our calculation is carried out for a quantum dot con- 
nected to two electronic reservoirs, denoted L and R, by 
single-channel leads. The reservoirs are characterized by 
their chemical potentials, [II and respectively, such 
that the electronic population in each of them is given 
by the Fermi distribution 



fa(E) 



1 



= (£-^)AbT + i 



(2) 



Below, we measure energies relative to the common Fermi 
energy, 



£ f = +Vr) 



(3) 



The bias voltage V across the quantum dot is 

V = (Ml - M R )/e . (4) 

One may consider two types of current-current corre- 
lations. The correlation of the net current, 

i = (h- 4)/2 , (5) 

flowing through the system, C(~\ui), is [see Eq. (JTJ] 



(6) 



C [ -\u)= I dte lbJt (6I(t)5I(0)) 

C (auto) H -C (x) (u;) 
Here we have introduced the definitions 

C^°Xw) = \(c LL {w) + C RR {u) 



(7) 



for the auto correlation and the cross correlation, respec- 
tively. [Note that both C^ nto \uj) and C<»(u;) are real.] 
The other correlation function is associated with the rate 
by which charge is accumulated on the dot, and is given 
by 



C {+) (uj)= j dte luJt (AI(t)Ai(0)) 

J —oo 

= i(c (auto) H + C (x) (w)) , (8) 



B. The noise spectra in the scattering formalism 

As discussed in Sec. UJ we treat the electronic inter- 
actions by the self-consistent Hartree-Fock approxima- 
tion. In this approximation the Hamiltonian becomes 
effectively a single-particle one, and therefore one may 
use the (single-particle) scattering formalism ™^ to ex- 
press the current correlation functions, Eqs. (JSJ)-©, in 
terms of the elements of the scattering matrix. 

In the scattering formalism the system is prresented 
by its scattering matrix; when the system is coupled to 
two reservoirs by single-channel leads the latter is a 2 x 2 
matrix, denoted S aa >(E), which is a function of the en- 
ergy E of the scattered electron. One finds^ (we employ 
units in which h = 1 hereafter) 



C aa , (w) 



77' 



x/ 7 (£) l-/ y (£ • 



(10) 



where / is the Fermi function, Eq. $2$, and 



77' 



(E,u } )=A^,(a,E,E + u } )A 1 , 1 (a , ,E + u J ,E) , 

(11) 



with 



A yY (a,E,E') 



S*(E)S ay (E') . (12) 



The various correlation functions are each given as a sum 
of four processes: the intra-lead contributions, which 
include the combinations /l(_E)[1 — fh{E + u)} and 
f R (E)[l — fn(E + uj)], and the inter-lead ones, which 
contain f L {E)[l - f R (E + u)} and f L (E)[l - f a (E + u)]. 
The relative contribution of each process to C^\uj) is 
determined by the respective product of the Fermi func- 
tions. In particular, at T = 0, this product vanishes 
everywhere except in a finite segment of the energy axis, 
determined by the bias voltage [see Eq. (Q}] and the fre- 
quency lo. Finite temperatures broaden and smear the 
limits of this segment, while varying the bias voltage may 
shift it along the energy axis or change its length. 



C. The self-consistent Hartree-Fock approximation 

Here the model Hamiltonian we use is specified, and 
the self-consistent Hartree-Fock approximation by which 
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it is reduced to an effectively single-particle one is de- 
scribed. That effective Hamiltonian is employed to derive 
the scattering matrix. 

Our model system consists of a multi-level quantum 
dot coupled to two (spinless) electronic reservoirs via two 
single-channel leads. The Hamiltonian of such a system 
is 

^ = W-leads + ^-tun + ^-dot i (13) 

where the leads are described by 

U-leads = E E £ ka c ka C ka 7 (14) 
a=L,R k 

with £fe Q being the energy of an electron in lead a, and c\ a 
(c ka ) is a creation (annihilation) operator of an electron 
with momentum k in lead a. The tunneling part of the 
Hamiltonian is 



E E 

\—L,R k,n 



(15) 



where Vk an is the tunneling probability amplitude of an 
electron to hop from lead a to the nth level on the quan- 
tum dot [<v n (d n ) is a creation (annihilation) operator of 
an electron on that levell . The Hamiltonian of the dot is 



H. 



dot 



E< 



: d) d 



u 



where e„ is the single-particle energy on the nth level of 
the dot. The electronic interactions (taking place solely 
on the dot) are described by the last term of Eq. (p~6|) . 
adopting the "orthodox" form^i in which U = e 2 /C rep- 
resents the capacitive charging energy (C is associated 
with the capacitance of the dot) and J\f g is proportional 
to the gate voltage. The role played by this parameter is 
discussed below. 

As is mentioned above, the interactions are treated 
within the framework of the self-consistent Hartree-Fock 
approximation^ This is accomplished as follows. The 
product of four operators appearing in Eq. (|16p . assum- 
ing n ^ m, is replaced by 



— {d n d m )d\ n d n — d n d m (d m d n ) 

-(did n )(did m ) + \(did m )\ 2 



(17) 



The decomposition given by the first two terms on 
the right-hand side of Eq. (IT7|) represents the Hartree 
approximation^ (these are also referred to as "diago- 
nal terms" ) while the next two terms represent the Fock 
one (and are also referred to as "non-diagonal terms"). 
The Fock terms were discarded in some of the previous 
works^S though in Ref. [3(3 they were shown to be of 
paramount importance. We show below that this is not 
the case in our calculation. 



In the Hartree-Fock approximation, the dot Hamilto- 
nian (now denoted Hdot) becomes effectively a single- 
particle one, 



U 



do t 



n n,m 



(18) 



where constant terms have been omitted. In this form of 
the Hamiltonian, the single-particle energy, e n , is 



while the effective hopping amplitudes, J nm , are 

Jnm = U(d m d n ) . 



(19) 



(20) 



These new parameters of the Hamiltonian are calculated 
self-consistently. As our system possesses time-reversal 
symmetry, the effective hopping amplitudes J nm can be 
assumed to be real (see below). 

For the sake of clarity, we confine ourselves to very 
low temperatures and to the linear-response regime, and 
consequently carry out the self-consistent calculation of 
~e n and J nm at zero bias, using the relation 



Am) o_ 
Z7T 



dEf(E)[g r mn (E) - (g r nm {E) 



The retarded Green function matrix of the dot, 

-l 



G 7 nm(E) = 



E~H 



dot 



T, r (E) 



(21) 



(22) 



is expressed in terms of the effective Hamiltonian (|18j). 
and the self-energy matrix E r due to the coupling with 
the leads is given by 



cm Kam 



E - e ka+ i7 l 



(23) 



with ?7 -> 0+. In Eq. (|2TJ f{E) denotes the thermal- 
equilibrium Fermi function of the leads obtained from 
Eq. © by using u a = ep. As usual, — lraI7 nn (E) de- 
scribes the broadening of the energy levels on the dot 
due to the coupling with the leads. 



D. The scattering matrix and the self-consistent 
calculation 

In order to obtain an explicit expression for the scat- 
tering matrix, it is convenient to present the tunneling 
Hamiltonian [sec Eq. (|15|) ] in the site representation, 
modeling each of the leads by a one-dimensional tight- 
binding chain of zero on-site energies and equal nearest- 
neighbor hopping amplitudes —j (for simplicity, we as- 
sume the two leads to be identical) . The coupling of the 
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left (right) lead to the nth level on the dot is denoted 
-Vhn (—Van)- One then has VkLn = —^/2/NV Ln smk 
and VkR n = — y2 /NVn. n sin k , where each of the two 
tight-binding chains consists of N sites whose lattice 
spacing is taken as unity. It follows that the scattering 
matrix at energy E = — 2j cos k is 



S(E) 



HE) \ v L g r {E)vl v l q t (e)vX 
f [v R Q r {E)vi v R g r (E)vl 



(24) 



where Vl (Vr) is the row vector of coupling amplitudes of 
the dot levels to the left (right) lead consists of vector el- 
ement V^/nn, and v(E) = 2j sinfc is the group velocity of 
the scattered electron. Time-reversal symmetry dictates 
a symmetric scattering matrix. This in turn implies that 
the Green function is symmetric in the indices n and m, 
and consequently (d^d^ [see Eq. (|2"Tj) ] must be real. 

In the following, we ignore the real part of the self 
energy of the dot Green function [see Eqs. (|2"2"|) and (|2"5|) ]. 
This is a plausible approximation provided that the real 
part does not vary significantly with the energy. Then, 
the (imaginary part of the) self energy consists of the 
coupling of the dot levels to the left and to the right 
lead, and reads 



12 r (E) 



T L (E) + T R (E) 



(25) 



where the (symmetric) matrices ^l(r) are given by 
v(E) 



L(R) 



(E) 



,■2 V L(R) 



V, 



L(R) 



(26) 



Employing these expressions in Eq. (|2"Tj) leads to 



00 dE 

—f{E) g r (E)ii; r (E)g a (E) 



, (27) 



with g® m (E) = g^* m {E). One may estimate the magni- 
tudes of the averages (e£jjd m ) by considering Eq. (|27]l in 
the limit of weak couplings to the leads, T/j <C 1, where 
f ~ vV 2 /j 2 denotes the typical value of the broadening 
of the dot levels due to that coupling. It is then expe- 
dient to use the Lehmann representation for the Green 
function, 



g a (e) ~ ( 



E - H dot ± if] 



E 



E±ir]-E e 



(28) 



where the Eg 's are the eigenvalues of Hd and the } 's are 
the corresponding eigenvectors. Notably, the Eis differ 
from the original onsite energies on the dot, e n , because 
of the Hartree corrections [see Eq. (fl"9)) ] and due to the 
Hartree-Fock ones [see Eq. (|20])]. which render the dot 
Hamiltonian non diagonal. However, when the coupling 
to the leads is small, the averages {d\- L d m ) for n ^ m 



are small as well, bringing Ei close to one of the original 
onsite energies and the eigenvector \ipg) to be close to the 
corresponding original basis vector for the noninteracting 
system. Inserting Eq. (|2"5)) into Eq. (|2T|) then yields 



E 



i 



dEf(E) 



\^i x )^h\i^L{E)+T R {E))\^){^ 
(E — + iij) (E — E„ — irj) 



(29) 



Examination of this expression (in the zero-temperature 
limit) reveals that in most cases the result of this integral 
is rather small, of the order of f/|e„ — e m \. A significant 
result is obtained only when n = m (and consequently 

These observations are borne out by numerically solv- 
ing the self-consistent equations (|2T|) . The computations 
have been carried out for a dot including 61 equally- 
spaced levels, e„ = nA (n = —30, ...,30), which are all 
coupled identically to each of the two leads^ That cou- 
pling [see Eq. (Ti"5j) ] has been modeled by a saddle-like 
point contact-potential^ 



r a (e F ) 



(l-^-T a (e F )) 2 , 



(30) 



where A is the level spacing, and a = L or R. The 
saddle-like point contact is characterized by its transmis- 
sion probability^! 



T a (e F ) = 



1 



1 



»)/" 



(31) 



where uj a and e a are the parameters characterizing each 
of the point contacts. 

Representative results of the self-consistent computa- 
tions are portrayed in Figs. [1] and [21 as a function of the 
interaction U. To produce those curves, we have chosen 
parameters such that the energy levels on the dot remain 
well-separated also after being modified by the Hartree 
terms and the Fermi energy is chosen such that for U = 
the energy levels are either occupied or unoccupied. In 
the lower panel of Fig [T] we see that lines representing 
the effective energy levels, Eq. (fT9")) . start linearly, since 
(d^dn) are either zero or one, as shown in the top panel 
of Fig. [T] As soon as a level crosses the Fermi level (the 
n = level), its occupation change from 1 to (see the 
top panel) and as a consequence the linear slope of the 
effective energies changes. After the second level (the 
n = 1 level) crosses the Fermi level the levels below the 
Fermi level return the their non-interacting values and 
the levels above the Fermi level grow linearly with U. 
The parameter N g was chosen to be 30 in FigQ] We find 
that in this parameter range, the Fock averages are neg- 
ligible almost always. It is only when the nth modified 
level, crosses the Fermi energy (upon varying U) that 
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FIG. 1: (color online) The top panel depicts the averages, 
computed self-consistently from Eq. (|21|). as a function of 
the interaction strength U. The horizontal curves are the di- 
agonal averages, (d^d„); the "cloud" around the lowest ones 
represents all the off-diagonal averages (d£d m ) with n ^ m. 
The red (thin-dashed) step-like curves represent the diagonal 
average corresponding to the energy levels crossing the Fermi 
energy ep , computed within the Hartree-Fock approximation; 
the purple (thick-dotted) curves show the same average com- 
puted within the Hartree approximation alone. Lower panel: 
the effective energies e n , Eq. (|19[) . are presented by the blue 
(dashed) and red (dashed-dotted) curves. The green (contin- 
uous) curves are J nm , Eq. (|20|) . and the dotted (black) hori- 
zontal curve marks the Fermi energy. Parameters: 6f = 0.5 ; 
<ll = £r = 0.598 ; uj l = lor = 0.8 ; M B = 30. The level widths 
[see Eqs. (30) and (3D] are then equal, T l =Tr = 0.05. (All 
energies are in units of the level spacing A.) 

the non-diagonal averages {(t n+i d n+i ,) , of small integer 
values i and i', acquire significant values. A further in- 
vestigation of the Fock terms is given in the Appendix 
for a two-level dot. 

The only parameter distinguishing Fig. [2] from Fig. [T] 
is Af g , taken to be 33 in the former. The parameter Af g 
determines the number of levels below the Fermi energy 
at large interaction. In Fig. [5] we see a similar behavior 
to the one in Fig. Q]only that now the effective energy 
levels have a negative slope for small values of U (as 



FIG. 2: (color online) Same as Fig. [TJ with j\f g = 33. 



shown in the lower panel) and that the total number of 
occupied levels in the dot increases as a function of U 
(as shown in the top panel). Denoting the total number 
of energy levels on the dot by 2nd + 1, we find that for 
Af g — iid > ep/ A [Af g — rid < £f/ A] , the modified energy 
levels shift as a function of U such that at large J7's there 
will be mod(Af g ) [mod(Af g ) + 1] levels below the Fermi 
energy. This is indeed the case in Figs. Q]and[2l 



III. THE NOISE SPECTRUM 

As the self-consistent determination of the effective dot 
Hamiltonian, Eq. (TT8l) . is strictly-speaking valid only at 
zero bias and zero temperature, we present in this section 
the noise spectra calculated in this regime. We comment 
on the effects caused by a finite bias voltage and finite 
(but low) temperatures in Sec. IIVI Our computations 
are confined to the case in which all levels on the dot are 
coupled identically to each of the two leads. In that case, 
the scattering matrix takes the simple form [see Eqs. 



7 



and (26])] 

S(E) = -1 

+ ig r (E) 
where 



T L (E) T^ 2 {E)Y R /2 {E) 



g r {E) = Y,Gnm(.E) 



(32) 



(33) 



depends on the sum of the two widths [see Eq. (|25|) ] 

T{E)=Y L {E)+T R {E) . (34) 

Furthermore, since the computations below are confined 
to the linear-response regime, Tl(E) and Tr(E) are re- 
placed by their values at the Fermi energy, as given in 
Eqs. ((301) and (H]) • It therefore follows that 



g r (E)-g a (E) = -iTg r (E)g a (E) 



(35) 



Using Eqs. (|32|) and (|35l) in conjunction with Eqs. ([7]) 
yields 



p 2 f°° 

C (auto) (w) / dEfWll-fiE- 

J — oc 

x (r 2 [g r (E)\ 2 + \g r (E + u;)\ 2 } 

-[rl + r^^^^ + ^ + cc] 



(36) 



and 



e 

2^ 



d£7/(£0[l-/(^ + w)] 



xr £ r fl ( s r CEOff°(£ + w) + cc 



(37) 



from which the two correlations introduced in Sec. Ill BL 
C ( -^ [Eq. © for the net current] and C (+) [Eq. ©, 
for the charge fluctuations] can be derived. Interestingly 
enough, the latter correlation depends only on the sum 
of the two widths, T, and therefore is insensitive to the 
(possible) asymmetry of the couplings of the dot to the 
leads. 

We begin our analysis with a discussion of the zero- 
temperature dc conductance, as given by the Landaucr 
formula 



G 



2tt 



T(e F ) 



(38) 



with the transmission coefficient T(ef) = rLr_R|<? r (eF)| 2 . 
In the absence of the interaction, the conductance as a 
function of the Fermi energy exhibits peaks whenever the 
Fermi energy coincides with one of the levels of the dot, 
see the solid line in Fig. [31 The peaks of the conductance 
broaden as the Fermi energy increases; this results from 
the energy dependence of the coupling between the leads 



and the dot, modeled by point contacts [see Eqs. ((30)) 
and (f3"Tj) ]. When the interactions on the dot are taken 
into account (in the self-consistent HF approximation) 
the levels on the dots are shifted and so are the conduc- 
tance peaks (see the dashed curve in Fig. [3]). Moreover, 
when the transmission probability of each of the point 
contacts is small, the dot is in the so-called Coulomb 
Blockade regime. This blockade can be lifted by tuning 
the Fermi energy or the gate voltage (equivalent to the 
parameter Af g in our model). In our case, the lifting oc- 
curs when the Fermi energy of the electrons in the leads 
coincides with one of the onsite levels on the dot; this 
happens at energy intervals of size A + U, and is associ- 
ated with the addition of one electron to the dot. 




FIG. 3: (color online) The dc conductance as a function of the 
Fermi energy. The continuous (blue) curve is the conductance 
of a non interacting dot, Eq. I)38[l; the dashed (red) curve 
corresponds to U — 0.5. Parameters used in Eq. pip : ujl = 
lor — 0.8 ; €l — — 0. (All energies are in units of the level 
spacing A.) 

Turning now to the noise spectra, we portray below 
our results for the net current correlation, Eq. ©. Fig- 
ure [4] shows the noise spectra (as a function of the fre- 
quency lu) for various values of the interaction strength, 
for a dot coupled symmetrically to both leads. The left- 
most curve in the main panel depicts the correlation in 
the interaction-free case. The spectrum shows a series 
of steps, of height V (in units of e 2 ), located roughly 
at frequencies corresponding to the renormalized onsite 
energies on the dot. As the interaction strength is in- 
creased the levels are shifted due to the Hartree correc- 
tion. (Those effective onsite energies are depicted in the 
inset of Fig. 01) As a result, the spectra are shifted as 
well. This step structure is similar to the ones obtained 
before, in the case of a non interacting two-leveli 2 - and 
sine le-level 13 dots. The curves in Fig. [4] are calculated 
within the Hartree approximation alone; including the 
Fock terms almost does not change them (see below). 

When the coupling of the dot to the leads is not sym- 
metric, i.e., Tl 7^ Tr, there also appear dips in-between 
the steps, as seen in Fig. [5] The frequencies at which the 
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FIG. 4: (color online) The noise spectrum for a dot coupled 
symmetrically to the leads, as a function of the frequency, for 
various values of the interaction strength (U = 0, 0.1,... , 
1, from left to right) computed within the Hartree approxi- 
mation. Inset: the effective single-particle energy levels £„/ A 
on the dot, as a function of the interaction strength for the 
same parameters. Parameters: €f — —0.8, 7V 9 = 30 and 
T L = T R = 0.025. These widths are given by Eqs. (J30j and 
(|31[1 with e_L = £fl = and uil = (jJr = 0.8. (All energies are 
in units of the level spacing A.) 



dips are located correspond to the spacing between any 
two levels that happen to be on the two sides of the Fermi 
energy, as seen in the inset of Fig. [5] where the effective 
energies are plotted as a function of the interaction. We 
discuss further the origin of the dip structure in Sec. IIVI 
As mentioned above, we have chosen our parame- 
ters such that the onsite levels on the dot remain well- 
separated also when modified by the interaction; as a 
result of this choice, the Fock terms' contribution is al- 
most negligible, and therefore the results shown in Figs. 
S]and[5]were computed within the Hartree approximation 
alone. To further clarify this point, we compare in Fig. [B] 
the correlation C^~\lo) obtained in the interaction-free 
case with the one computed by the Hartree approxima- 
tion, and within the Hartree- Fock one. As can be seen in 
that figure, the last two curves are almost indistinguish- 
able. 

Notwithstanding the significance of studying the noise 
spectra as a function of the interaction strength, it seems 
more practical to explore it when the Fermi energy (or the 
lead chemical potential at finite, but low, temperatures) 
is varied. This dependence is presented in Fig. [7J for an 
interaction-free dot. The dependence of the spectrum on 
the Fermi energy originates from the ep — variation of the 
point contact's transmission. As the latter is decreasing, 
the height of the step seen in the spectrum is decreasing. 
The location of the step in Fig. [7J corresponds to the dif- 




FIG. 5: (color online) Same as Fig. [4]for a dot coupled asym- 
metrically to the leads. The vertical red line in the inset indi- 
cates the chosen interaction strength, (7 = 0.5, and the dashed 
black line indicates the Fermi energy, ef = —0.9. The dot lead 
coupling parameters are: Vt = 0.022 and Fr — 0.014. The 
widths are found from Eqs. (|30[) and (|31[) with ei, = 0, €r = 1 
and ujl = uir = 0.8. 



0.4 




co/A 

FIG. 6: (color online) Same as Fig. [4] for an interaction- 
free dot, shown by the solid (black) curve, and for U = 0.7 
computed within the Hartree approximation, shown by the 
dotted-dashed (blue) curve, and within the Hartree-Fock ap- 
proximation, shown by the dashed (red) curve. The Fermi 
energy is e.p = —0.7. 



ference between the Fermi energy and the closest energy 
level. We have chosen parameters for the Fermi energy 
in Fig. [TJsuch that the different steps will appear around 
the same frequency. 
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FIG. 7: (color online) The noise spectrum of a non interacting 
dot as a function of the frequency, for various values of the 
Fermi energy, ef = —0.7, —1.3, —1.7, and —2.3 (from left to 
right). Parameters: r L = F R = 0.028, 0.014, 0.009, 0.004 
[from left to right, obtained from Eqs. (|30[1 and (|31[) with 
ojl = = 0.8 and £l = £r — 0]. (All energies are in units 
of the level spacing A.) 

IV. DISCUSSION 

We have presented a calculation of the noise spectra 
of the currents passing through a multi-level dot, treat- 
ing the Coulomb interactions (assumed to take place on 
the dot alone) within the self-consistent Hartree-Fock ap- 
proximation. The computations are carried out for the 



noise associated with the net-current fluctuations, as this 
quantity is intimately related to the ac conductance. Our 
main conclusion is that this correlation, as a function of 
the frequency, reflects the energy spectrum of the dot. 
It shows steps at frequencies corresponding to the onsite 
energy levels of the dot, and dips, occurring at frequen- 
cies corresponding to the spacings between any two levels 
located at the opposite sides of the Fermi energy. The 
latter appear only when the dot is coupled to the leads in 
an asymmetric form. This pattern is similar to what has 
been obtained befor o 12 ' 13 in the case of a few-level non 
interacting dots; this is not very surprising since we have 
operated in the parameter regime where the Coulomb in- 
teraction did not modify drastically the energy levels of 
the dot. For technical reasons, the computations were 
restricted to zero temperature and zero bias voltage. 

It is of course very desirable to examine the effects of 
finite bias voltages and finite (though rather low) tem- 
peratures. In this section we carry out such a study, con- 
fining ourselves to the simplest case of a two-level non 
interacting dot. In the process, we also gain more in- 
sight into the origin of the step and dip pattern of the 
noise spectra. In particular we focus on the dip structure 
since (as mentioned in Sec. fl]) it is debated in the liter- 
ature. We believe that the conclusions drawn from this 
simplistic model are valid, at least qualitatively, also for 
the interacting multi-level dot, as long as the the inter- 
action is not too strong and the number of onsite levels 
is sufficiently large. 

We begin with the correlations C^> at zero bias and 
zero temperature. Then, Eqs. (136[) and (|37[) in conjunc- 
tion with Eqs. © and © yield 



-in 



dE^T, 



(\9 r (E)\ 2 + \9 r (E + ^)| 2 ) - (T L ± r R f(g r (E)g"(E + u ) + cc) } 



(39) 



r 



In the case of a two- level, non interacting dot, the Green 
function Eq. f|33|) takes the form 



9 r (E) = 



(E - e,) + (E - e 2 ) 



(E-e 1 +i*k)(E-e 2 +i^) + 



r- 



(40) 



In particular, when the two levels are well-separated, 
g r consists roughly of two Breit-Wigner resonances of 
(equal) widths T/2, centered around ei and £2. The con- 
tribution of the first term in Eq. (|39[) to the integration 
is hence quite clear: As long as none of the resonances 
is within the integration bounds the integral will almost 
vanish, while whenever the frequency integration encom- 
passes one of the Breit-Wigner resonances, due to g r (E) 
[g r {E + u)] in case the energy level is below [above] the 
Fermi energy, the result will be a step in the noise curve, 
whose height is roughly T/2 (in units of e 2 ). The second 



term in the integrand of C*- - ) alone is non zero only when 
the dot is coupled asymmetrically to the leads, namely, 
when Tl ^ Tr. (This is not the case for the charge- 
related correlation C^ + K) The integrand of this term 
consists of products of resonances related to energy lev- 
els below and above the Fermi energy; the contribution to 
the integration will be significant whenever the frequency, 
oj, will roughly match the level spacing, e\ — £2 — w. This 
contribution is negative, and leads to the dip in the cor- 
relation as a function of the frequency of depth roughly 
equal to (Tl — Tr) 2 /2T (in units of e 2 ). Interestingly, 
the dip in suppresses the noise almost completely. 

Figure[8]depicts the behavior described above. In both 
correlations, and C^~\ there appears a dip at the 

frequency corresponding to the difference ei — £2; how- 
ever, in the case of C(~\ that dip disappears once the 
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FIG. 8: (color online) The correlations C' + ' (upper panel) 
and C' - ' (lower panel) for the two-level case, as a function of 
the frequency at zero temperature and zero bias. The different 
curves in the lower panel correspond to Tl — 0.1 and Tr = 
(brown continuous curve), Tl = 0.0745 and Fr = 0.0255 
(blue dashed curve), and Tl = Tr = 0.05 (black dotted 
curve). Parameters: ei = 1, e 2 = —2, and ef = 0.5. All 
energies are in units of A. 



dot is coupled symmetrically to the leads. In contrast, 
the correlation C"~) is not sensitive to this asymmetry, 
as expected. The behavior portrayed in Fie. [5] is in 
agreement with the findings of Refs. HH andQ, but con- 
tradicts those of Ref. [H, where a dip was obtained even 
in the case where the dot is coupled symmetrically to the 
leads. The two levels discussed in Ref. are also cou- 
pled to one another, making it a different system than 
ours. The Green function g r [c.f. Eq. describing 
that geometry is 



g r (E) = 



{E-e 1 ) + (E-e 2 )-2J 
(E - ei + i$)(E -e 2 + il) - (J + iT) 2 



(41) 



where — J is the coupling between the levels. Since the 
form of Eq. is unaltered by the levels coupling [one 
has just to employ Eq. f4"Tj) ] a dip will not appear in the 
noise for Tl = Tr. 

At finite biases Eq. is not applicable anymore and 
one needs to use Eqs. (fTUjl - p^]) to obtain the different 
contributions of the inter- and intra-lead processes to the 
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FIG. 9: (color online) The noise spectra C' + '(o;) (top figures) 
and C^~'(oj) (bottom figures) of the biased dot at zero tem- 
perature with ei = A and £2 = — 2A. The Fermi energy is 
set to 0.5 in all figures. The different figures are for for eV=3 
(first and third figures from the top) and 6 (second and forth 
figures from the top). The five curves correspond to Tl —0.1 
and r_R = (brown continuous curves), Tl = 0.0745 and 
T R = 0.0255 (blue dashed curves), T L = T R = 0.05 (black 
dotted curves), T L = 0.0255 and T R = 0.0745 (gray long 
dashed curves), and Tl = and Tr = 0.1 (green dotted- 
dashed curves'). All energies are in units of A. 
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noise spectra. As pointed out in Sec. Ill B[ a finite bias 
modifies the integration bounds in Eqs. (fTUl) - ([T^t . and 
as a result the steps in the noise spectra are not only 
shifted, but split as well, as can be seen in Fig. [§J The 
dip however, (which occurs roughly at oj = |ei — ea|) is not 
affected by the finite bias, in agreement with the findings 
of Ref. [l6j. When only one of the levels lies in-between 
the two chemical potentials of the leads (this scenario 
corresponds to eV — 3 for the parameters used in Fig. 
[5]) the dip in becomes sensitive to the ratio F^/Tr, 
and the dip in C*- - -* can turn into a peak for certain values 
ofTi/r^. When the bias is large enough, so that the two 
levels are well between the two chemical potentials of the 
leads (for the parameters used in Fig. [9j this occurs for 
eV = 6), the dip in no longer depends on the ratio 
Tl/Fr, and the dip in G^' turns into a peak. Finally, 
when the bias is large enough for the noise to become non 
zero at a; = — |ei — £2!, an additional dip [peak] appears 
in C^ + ' [C- - - 1 ] at negative frequencies. 

The effect of both a finite bias voltage and a finite tem- 
perature is exemplified in Fig. 1101 Interestingly enough, 
while the steps in the spectra are smeared by the tem- 
perature, the structure around (J = \e\ — £2! seems to be 
rather immune. The dependence of this (dip) structure 
on the asymmetry of the dot couplings to the leads may 
be therefore amenable to experimental verification, mak- 
ing the noise spectra a possible tool to study the energy 
spectrum of quantum dots. 



C (+ \(o)/e 2 




FIG. 10: Same as Fig. H with T = 1.5A and eV = 6A. 
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particular its dependence on the off-diagonal coupling of 
the levels. By the Landauer formula,^ the differential 
conductance is 

J^F] = ^(^tld^l) + T tld (^)) , (Al) 



Appendix A: A possible experimental observation of 
the Fock terms 

As can be appreciated from Eq. (fT5|) . the Fock terms 
couple together the (Hartree-modified) single-electron 
energies on the dot. Thus, the original onsite levels are 
modified by this coupling. A possible way to detect these 
terms might be to monitor the outcome of such a cou- 
pling between the dot's levels. For example, one may 
ask how the I-V characteristics are modified when this 
off-diagonal coupling is changed. This may be related 
to the recent experiment of Konemann et al.^ in which 
the slope of the steps in the differential conductance de- 
pended on the polarity of the bias voltage. These authors 
explain their results by relating them to the broadening 
of the levels (see below). 

In order to explore this possibility, we consider the sim- 
plest case of a (non interacting) two-level dot and calcu- 
late its differential conductance at zero temperature, in 



dI/d(eV) 




FIG. 11: (color online) The peak in the differential conduc- 
tance as a function of the bias voltage, for various couplings J: 
J = (black continuous curve), J = 0.1 (red dotted curve), 
and J = 0.2 (blue dashed curve). Parameters: e-p = 0.5, 
ei = 1, 62 = 2, Tl ~ r_R = 0.05. (All energies are in units of 
the level spacing A.) 
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where 7tld is the transmission of the two-level dot. De- 
noting the two (original) onsite energies by ei and £2, 
their coupling by —J, and assuming that they are iden- 
tically coupled to each of the two leads, the transmission 
is 

T ThD (E)=T L T R \g r (E)\ 2 , (A2) 

where r = T L +T R and g r (E) is given in Eq. flU). (Note 
that we have assumed here that the widths T l and Yr are 
independent of the energy.) The differential conductance 
is plotted in Fig. [TTJ As can be seen there, increasing 
the coupling J between the levels narrows down the peak 
of the differential conductance, in a rough agreement 
with the observations and interpretations of Kdnemann 



et al££ Indeed, as can be noted from Eqs. fill) and (|A2I) . 
the off-diagonal coupling J modifies the broadening of 
the energy levels. In the simple case treated here, the 
width of the upper level (say ei) is reduced, from T to 
r(l — b), while the width of the lower level is enhanced, 
to be r(l + b), with b = 2J/^AJ 2 + (e 1 -e 2 ) 2 . This ob- 
servation helps us to understand the (small) differences 
between the curves in Figs. [T] and [2] calculated within 
the Hartree approximation alone, and those computed 
within the Hartree-Fock approximation. Consider for ex- 
ample the diagonal averages, corresponding to the level 
occupancies. As is known,— their variation with the in- 
teraction strength is determined by the levels' width; the 
Fock terms induce modifications in those widths and con- 
sequently change slightly the curves. 
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